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Abstract 

Equations related to the Polchinski version of the exact renormalisation group equa- 
tions for scalar fields which extend the local potential approximation to first order in a 
derivative expansion, and which maintain reparameterisation invariance, arc postulated. 
Reparameterisation invariance ensures that the equations determine the anomalous di- 
mension T] unambiguously and the equations are such that the result is exact to O(e^) in 
an e-cxpansion for any multi-critical fixed point. It is also straightforward to determine 
7] numerically. When the dimension d = 3 numerical results for a wide range of critical 
exponents are obtained in theories with 0{N) symmetry, for various N and for a ranges 
of rj, arc obtained within the local potential approximation. The associated rj, which 
follow from the derivative approximation described here, are found for various N. The 
large N limit of the equations is also analysed. A corresponding discussion is also given 
in a pcrturbative RG framework and scaling dimensions for derivative operators are 
calculated to first order in e. 
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1 Introduction 



Exact functional renormalisation group flow equations [1, 2, 3, 4, 5, 6, 7, 8] allow, at least 
for scalar field theories, the possibility of a non perturbative analysis of fixed points and 
determination of critical exponents which control the RG fiow near any fixed point. In 
all such equations there is a cut off^ function K{p^) which is essentially arbitrary save for 
i^(0) = 1 and vanishing sufficiently rapidly as oo. Any physical results, such as 

precise values for exponents, should be independent of the cut off although it may be 
feasible to optimise over different cut off functions [9]. The exact RG flow equations are 
hard to handle except in some truncation or expanding in perturbation theory. The local 
potential approximation (LPA) neglects the spatial dependence of the flelds ^ and reduces 
the effective action from a highly non trivial functional to a simple function V{<\)) 
and the RG flow equations become a non linear differential equation for V of the form 
Vt = F^t; yl-, yl') where V denotes the derivative with respect to t = — log A, with A the 
cut off scale (the equations are invariant under rescalings of A). The potential RG flow 
fixed points, Vt — F as t — oo, are determined by requiring smooth solutions for all (f) of 
F(y, V', V") = 0. The critical exponents describing the RG fiow in the neighbourhood of 
the fixed point, Vt = y + e'^' *?;,., may then be calculated by finding the eigenvalues 

Xr, and eigenfunctions Vr, for a corresponding linear differential operator depending on the 
fixed point solution for V. For all > the associated operators are relevant and it is 
necessary to tune = to attain the fixed point under RG fiow. 

When applied to the Polchinski RG equation [4], for which F{V,y' ,y") has a very 
simple quadratic form, the LPA has the virtue that all dependence of the cut off can be 
removed by rescalings of V and (j). Although rather crude the LPA is compatible with the 
global features of RG flow since, in cases that have been investigated, it realises the same 
fixed points as are present in the full quantum field theory for scalar fields that are found 
by other techniques (this is not manifestly true for more complicated theories with gauge 
fields and fermions [10]). 

Despite describing the essential features of the landscape of critical points for scalar 
theories the LPA has, nevertheless, many limitations. In particular it is not possible to 
consistently determine r], the anomalous dimension for (p. In theories with dimension d = 3, 
r] is generally small but clearly results for critical exponents must then have an error of 
at least 0{rj), although results when the LPA is applied to different RG flow equations 
differ in general by rather more than this. Attempts to go beyond the LPA usually invoke 
an expansion in terms of derivatives of (p [11, 12, 13, 14, 15, 16, 17, 18]. To flrst order 
this introduces a function Z{(j)) which is the coefficient of {dcf))"^ in an expansion of the 
effective action (for multi-component fields (pi this becomes a symmetric tensor Zij{(p)). Z 
and y obey coupled equations which in principle allow rj to be determined by requiring 
non singular solutions for both Z and y . However the dependence on the cut off becomes 
more severe in the derivative expansion. Applied to the Polchinski equation there are two 
constants A, B which are essentially arbitrary [13]. Apart from this arbitrariness the results 
also depend on the value chosen for Z{Q) in solving the coupled equations [15]. 

Exact RG equations, without approximations, are invariant under reparameterisations, 
including rescalings, of the fundamental fields [19]. This property ensures that the full 
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equations have a line of physically equivalent fixed points which may be parameterised by 
different values of Z{0), [15, 20, 5]. Physical results, such as rj, are independent of where on 
this line the fixed point solution is chosen. As a consequence of the line of equivalent fixed 
points the calculated exponents must include one which is exactly zero. The corresponding 
marginal operator is redundant, essentially one which vanishes on the equations of motion. 

In the context of the Polchinski RG equation it was shown, for arbitrary dimensions d, 
in [21] that for any local operator O, such that f d'^x O gives an eigen-operator represented 
by an eigenfunction for the linearised equations with critical exponent Xq , then it is possible 
to construct associated redundant operators with exponents 

A = Ao-i(d + 2-r/)-2m, m = 0,l,2, (1.1) 

irrespective of any particular choice of a smooth cut off function. Furthermore the operator 
is a local operator determining an eigenfunction with A^^ = ^(d + 2 — ry) and hence, 
applying (1.1) with m = 0, this directly shows that A = is a possible eigenvalue whose 
eigenfunction generates the marginal operator necessary for reparameterisation invariance. 

Although reparameterisation invariance is a property of the full non linear RG equations 
it is generally lost in approximations such as the derivative expansion. There is no longer a 
fluctuation eigenfunction with A = exactly. Here we heuristically construct equations for 
V {(/)), Z {(f)) , and also for multi-component generalisations, which maintain these desirable 
features. The equation for V remains the same as in the LPA except for the introduction of 
rj. The associated equation for Z depends on the solution for V and determines rj. Using an 
appropriate scalar product an integral expression may be found which may be used to find rj 
in a fashion which is manifestly independent of Z(0). The eigenvalues for the corresponding 
differential operator are in accord with (1.1) when m = and the zero mode eigenfunction 
can also be found explicitly. 

For the purposes of comparison we also discuss results for derivative operators of the 
form Gij{4>)d^(t)id^4>j using standard perturbation theory techniques and the e-expansion to 

obtain results for the anomalous dimensions at the fixed point to first order in e. For such 
derivative operators it is necessary to take account of mixing with scalar operators -F(0) with 
the same dimension but there are also additional constraints on the associated /3-functions. 
Keeping only contributions just to first order in the coupled {F, Gij) then an infinitesimal 
variation 5(j)i = Vi{(j)), for non-linear Vi{(/)), in the lagrangian is equivalent to corresponding 
changes in {F, Gij). This leads to identities which show that the scaling dimensions satisfy 
relations of the same form as in (1.1) in general in a perturbative framework. We are also 
able to determine the scale dimensions to 0(e) at each of the multi-critical fixed points for 
scalar theories when 0{N) symmetry is imposed. Although such operators are irrelevant 
as far as RG flows they arc of course of interest in determining the spectrum of operators 
and scale dimensions in the theory at its critical points. 

In this paper we describe in the next section results for the simplest case of a single 
component field which corresponds to the Ising model and has been much discussed pre- 
viously. For d = 3 the equations are solved numerically and the associated eigenvalues 

determined for various values of rj. The appropriate value of rj necessary for a non sin- 
gular solution of the Z-equation is also found. In section 3, we extend the discussion to 
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multi-component fields, imposing 0{N) symmetry so that simple equations, of similar form 
to those considered in section 2, are obtained. The eigenfunctions are then 0{N) tensors. 
The irreducible representations are given by symmetric traceless tensors of rank I and the 
corresponding eigenvalues depend also om I. Numerical results are then given for various 
N and /. In section 4 we show how these equations may be solved in an e-expansion recov- 
ering perturbative results at the various possible non-trivial multi-critical fixed points as 
the dimension d is reduced. In section 5 we consider perturbatively the usual /^-functions 
in a loop expansion, extending results obtained in the single component case in [21]. In 
section 6 these results are extended to derivative operators and mixing effects taken into 
account. In a conclusion we make some more general remarks concerning the status of 
the equations discussed in this paper. Although they have been motivated by requiring 
that they share general properties of the exact RG equations they serve to show how these 
may still be maintained in quite simple approximations. Various calculational details are 
relegated to four appendices. In appendix A we show how the equations can be solved for 
large N and a formula for rj to 0(A^~^) obtained which is quite close to the exact large N 
result. In appendix B we give some details of the perturbative results for /3-functions that 
are used in sections 5 and 6. In appendix C we give a general discussion using dimensional 
regularisation of the consequences of invariance of the regularised theory under variations 
5(l)i = — In the final appendix D we give some details of the nearest singularities 
that are found numerically when the solution of the local potential approximation for the 
Polchinski equation is extended to the complex plane. 



2 Equations for a Single Component Field 

It is simplest to consider first a single scalar field (p corresponding to the universality class 
for the Ising model. At a fixed point the equation for V{(p) is 

V"{(f)) -l{d-2 + rj) cf>V'{(l>) - V'{(f))^ + dV{(f)) = . (2.1) 

This is just the standard LPA for the Polchinski equation including the anomalous dimension 
rj. In general this is set to zero as there is no mechanism to determine this from (2.1). The 
two trivial solutions of (2.1) are V{(f)) = 0, for the Gaussian fixed point, and V{(j)) = 
J (2 — rj){<fP' — |), for the high temperature fixed point. Non trivial solutions even in (p, so 
that V'{G) = 0, which are non singular for all and 

. „ 2d 

V{(t)) \{2 - r]) (ff + A(j)'i+^-v for large 4> , (2.2) 

depend on a precise choice for V^(0) = k which then determines A. Such solutions appear 
whenever d is reduced below 2n/ (n — 1) for n = 2, 3, . . . [22, 23]. The critical exponents are 
then determined from the eigenvalue equation 

Vf{<p) = {\-d)f{c(>), (2.3) 

with the differential operator 

^'=^-|(<i-2 + ,)4-2n*)^. (2.4) 
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It is easy to see, using (2.1), that 



W((/)) = -i(d + 2-r/)F'(0), P(/) = -i(d-2 + 7?)(^-2F'(</)), (2.5) 

and hence we may construct two exact odd eigenfunctions 

U{cl>) = l{2-v)(t>-V'{4>) , h = ^(d+2-r?) , fR{(t>) = V'{cl>) , \r = Ud-2+n) , (2.6) 

/il(0) corresponds to a redundant operator with Xr given by (1.1) with O the identity 
operator, which corresponds to the solution of (2.3) /{(f)) = 1 with X = d. 

It is important for our later discussion to recognise that V in (2.4) is hermitian with 
respect to a scalar product defined by 

if, 9) = J d<^e-K'^-^+'')^^-^^(^) mg{<P) , (2.7) 

so that 

{f,Vg) = {Vf,g). (2.8) 

Extending the RG equations to Z{(p) we propose that, in conjunction with (2.1), the 
associated equation at the fixed point 

Z"{cP) - i(d - 2 + r?) cpZ'{<P) - 2V'{(^)Z'i^) - 2V" {<P) Z {cj,) 

= {V- 2V"{ct>))Z{<t>) = r/ - ^ V"{<t>f . (2.9) 

Together with (2.1) this satisfies reparamctcrisation invariance so that rj is independent of 
the particular initial 2^(0), unlike the case for other analogous derivative expansion equa- 
tions. To ensure non singular solutions for all requires only a special choice for 77. Asymp- 
totically, for large ^, solutions of (2.9) have the form 

m - 4^ - ^ + cc^-^ . (2.10) 

2{d + 2) 2-7] ^ ^ ' 

In general the value of the asymptotic constant C depends on Z{G). 

The proposed Z-equation (2.9) is similar to the derivative expansion result in [13]. It 
differs in that the coefficient of the V" Z term is 2, rather than 4, and that on the right hand 
side there is a definite coefficient rather than an essentially arbitrary cut off dependent 
constant (in terms of the equations in [13] we are taking for the cut dependent constant 
B the precise value g^)- In respect of these terms (2.9) is identical with an analogous 
equation obtained in [21] using an expansion in terms of scaling fields which is similar in 
spirit to the derivative expansion. In the scaling field approach the corresponding coefficient 
is determined precisely essentially by those divergencies in two point amplitudes which 
are universal, i.e. renormalisation scheme independent. If G{x) is the cut off dependent 
propagator, then we have for the following products 

^G(x)" ~ -nan6\x) , ^G{xf^-^ ~ -(2n - 1)6^ dH\x) , (2.11) 
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for 

2n 

d = dn = (2.12) 

n — 1 

such that a„, 5„ are constants, independent of the cut off function and depending only on 
the large x behaviour of G{x). According to the results obtained in [21] 

h}L = (2 13^ 

a„2 2{dn + 2) ■ ^ ■ ^ 

The coefficients in RG equations such as (2.9) should not depend on the particular critical 
point, here labelled by n, but may depend on the spatial dimension d. Applying the e- 
expansion to (2.9) with the particular coefficient ensures, as was shown in [21] and also 
subsequently here, that r] is correct to 0{e^) for all critical points n = 2,3, . . . Although 
results such as these for 77 were obtained from Wilsonian RG equations as soon as they were 
first proposed, and were shown to be independent of the detailed cut off function [24], they 
are also identical, of course, with results from standard Feynman graph techniques which 
arise directly from the coefficients of the universal logarithmic divergencies for particular two 
point Feynman graphs. These logarithmic divergencies are equivalent to (2.11). In a sense 
compatibility with the e-expansion may be regarded as an optimal choice for such constants 
as B. However, in the scaling field derivation described in [21] there is no free constant to 
determine and agreement with the e-expansion is not imposed but follows automatically. 

In addition (2.9) differs from corresponding equations in [13] and [21] by the absence 
of a 'qZ term. Removing such a contribution is essential to obtain subsequent results. In 
general in a derivative expansion there are also expected to be additional contributions on 
the right hand side of (2.1) involving Z but the exact form differs between [13] and [21] and 
also involves a cut off dependent constant which we arc here essentially setting to zero. 

Corresponding to (2.1) and (2.9) there are associated eigenvalue equations for critical 
exponents 

I ^ + d \(f^^A=x(f^^~^\ (2U) 

[-2Z'{<p)^-{2Z{<p)-^V"{cl>))^ V-2V"{ct>))\9{4>)) ' ^' ^ 

It is easy to see that this decouples into pairs of equations for eigenvalues \f,Xg where 
A/ is obtained from (2.3), with the corresponding g determined in terms of / by inverting 
V - 2V"{(I)) - Xf, and also, with / = 0, 

Vg{ct>)-2V"{4>)9{(I>) = >^9 9{^)- (2-15) 
For the eigenfunctions in (2.6) the corresponding functions g are given by 

-5^'^) = 9Ri^) = Z'icP) . (2.16) 

To show reparameterisation invariance of (2.9) we note that for any solution of (2.3) 
there is a corresponding solution of (2.15) given by 

= /'(</>), Xg = Xf-l{d + 2-r,), (2.17) 
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in accord with (1.1). Starting from f(f){(p) in (2.6) it is then easy to obtain an exact zero 
mode 

5o(0) = ^(2-r/)-y"(0), (2.18) 

representing the necessary marginal redundant operator present in the RG flow equations. 
Since T> — 2V"{(j)) is hermitian with respect to the scalar product in (2.7) we must have, for 
consistent solutions of (2.9), 

v{go,l)=V I d0e-i('^-2+'')^^-2^W(i(2-ry)-n'A)) 

d(/>e-i('^-2+'')^'-2^W (1(2 -7])- V"{c))))V"{(l)f . (2.19) 

Since r] is small it is easy to iterate (2.19) in conjunction with (2.1) starting from r/ = to 
determine the consistent solution for rj with high numerical precision. 

When d = 3 we may readily solve (2.1) numerically tuning k = V{0) so that the sin- 
gularity in the solution arises for the largest possible value of (p compatible with numerical 
precision, (2.1) was written as two coupled first order equations and were integrated from 
= using RK4. The limiting results when ry = are shown in Figure 1. Numerical re- 
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Figure 1: Typical numerical solution of (2.1), with d = 3,ry = 0, for T^(</>) starting from the 
critical value k = V(0) = 0.0761994008. Due to rounding errors the solution breaks down 
for (/> ^ 7 and is singular at (pQ ~ 11.23. As shown, for (f) S (po, it is well approximated near 
the singularity by the leading singular form V^ingX"^) = ~ ^T^iipo — (p) + \4>o{(po — (p) +const. for 
solutions of (2.1). Note that this has a minimum at (p = (pQ — 4/00 matching the minimum 
of the numerical solution. 

suits for the Z-equation (2.9) are also shown in Figure 2. These were obtained in a similar 
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fashion as for V in terms of corresponding first order differential equations. The solutions 
also develop singularities which are very sensitive to the value of rj, where the corresponding 
F-solution of course has been used. 



z 

1.6t 




(b) eta=0. 041347 

Figure 2: Numerical solutions of (2.9) for various Z(0) with rj just either side of the critical 
value so that the singularity arises for the largest possible 0. The graphs demonstrate how 
rj is independent of Z{0). 



Having determined V the eigenvalues are then determined numerically for small 
values of rj and c? = 3 by optimising the eigenvalue such that the eigenfunction blows up 
as slowly as possible withing the range where results for V{(/)) are reliable. The results are 
ordered such that An > An+i with Aq = 3. For n even, odd the associated eigenfunctions are 



7 



even, odd in (p. From (2.6) Ai = ^(5 — i]) and A3 = ^(1 + ??) wliich provides a consistency 
check on our numerical results. For even n the results are in Table 1 and for odd n in Table 
2. For the small values considered the dependence on r/ is close to linear. For 77 = our 
results agree with the much more accurate determinations in [25, 26]. 





A: 


A2 


A4 


Ae 


As 


Aio 


0.00 


0.076199401 


1.5395 


-0.6557 


-3.1800 


-5.9122 


-8.7961 


0.01 


0.073512228 


1.5409 


-0.6449 


-3.1560 


-5.8735 


-8.7416 


0.02 


0.070885912 


1.5421 


-0.6341 


-3.1319 


-5.8344 


-8.6866 


0.03 


0.068319137 


1.5433 


-0.6232 


-3.1076 


-5.7951 


-8.6312 


0.04 


0.065810847 


1.5444 


-0.6123 


-3.0832 


-5.7554 


-8.5754 


0.05 


0.063359963 


1.5454 


-0.6013 


-3.0585 


-5.7155 


-8.5190 


0.06 


0.060965439 


1.5463 


-0.5903 


-3.0337 


-5.6751 


-8.4621 


0.07 


0.058626258 


1.5471 


-0.5793 


-3.0087 


-5.6345 


-8.4048 


0.08 


0.056341435 


1.5478 


-0.5682 


-2.9835 


-5.5935 


-8.3469 


0.09 


0.054110012 


1.5485 


-0.5570 


-2.9582 


-5.5521 


-8.2884 


0.10 


0.051931059 


1.5490 


-0.5458 


-2.9326 


-5.5103 


-8.2294 



Table 1: Even Eigenvalues and initial value k = V{0) for non singular solutions 
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As 


A7 


Ag 


0.00 


-1.8867 


-4.5244 


-7.3377 


0.01 


-1.8696 


-4.4932 


-7.2911 


0.02 


-1.8524 


-4.4618 


-7.2422 


0.03 


-1.8351 


-4.4301 


-7.1970 


0.04 


-1.8177 


-4.3982 


-7.1493 


0.05 


-1.8002 


-4.3660 


-7.1012 


0.06 


-1.7826 


-4.3336 


-7.0527 


0.07 


-1.7648 


-4.3010 


-7.0038 


0.08 


-1.7470 


-4.2680 


-6.9545 


0.09 


-1.7290 


-4.2348 


-6.9047 


0.10 


-1.7109 


-4.2013 


-6.8544 



Table 2: Odd Eigenvalues 

We have also determined the value of rj required for non singular solutions of (2.9) and 

verified that the result is independent of the value chosen for Z(0) and also in agreement 
with (2.19). This may be used to determine r] by iteration starting from = and gives 

7/ = 0.041347 for A; = 0.0654776 . (2.20) 

For this value 

i/= I/A2 = 0.647, a; = -A4 = 0.612. (2.21) 
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3 Multi-Component Fields 



There are natural generalisations of the above equations to the case of a A'^-component 
scalar field ^j. Instead of (2.1) we have 

d'V{ct>) -Hid- 2)ct> + (r?</.)) • ^F(^) - • ^V{ct>) + dV{<i>) = , (3.1) 

for r]i.j a symmetric anomalous dimension matrix and = ^ ■ The corresponding 
equation for critical exponents is just as in (2.3), 

VF{<P) = {X-d)F{<l>), (3.2) 

but with the differential operator 

I> = #-l(«i-2)0 + (#))-|-2|m-|. (3.3) 

and we may now also allow tensorial eigenfunctions Fij...(^). Trivially PI = so that X = d 
is an exact eigenvalue. Just as in (2.6) we have exact vector eigenfunctions since 

VF^M = i(d - 2)F^44>) + i {rjF^i4>))i , F^44>) = l{n<j>)i - ^Vi4>) , (3.4) 

and we may choose a diagonal basis for r)ij. The corresponding scalar product to (2.7) is 
just 

{F, G) = j d^0 ^-\{d-2W~\Hn4')-2V{ci>) p^^^ ^^^^ ^ (3 5) 

with additional tensorial contractions if required. 
Extending (2.9) to Zij{(p) we require 

2d (9^ 

vz,,{4>) - 2V^,,{4>)z^),i<p) = r?,, - — v,kimd<i>) , y^Ji4>) = Q^.y(^) ' (3-6) 

with ( ) used to denote symmetrisation of the i,j indices. For fluctuations F,Gij around 
the fixed point solutions we have (3.2) and also the associated coupled equation 

VGiji<P)-2V^,ki<f>)Gj)kW 

- 2dkZ,j{cP) dkF{<t>) - (^2Z(,fc(</>) - J^y(^fc(0)j9.)5fcF((^) = AG,,(</>) . (3.7) 

As before it is easy to see that the possible eigenvalues are determined by (3.2) and \g 
given by 

VGij{<i>) - 2V^ik{<t>)Gj)k{<t^) = XcGiM) , (3.8) 

corresponding to (2.15). For any vector eigenfunction -Fj(^) there are associated eigenfunc- 
tions for (3.8) given in terms of 

Gij{<P) = diFji^) + djFii^) , (3.9) 
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since 

VGiji<t>) - 2V(,fc(<^)G,)fe(</.) = (A,. - lid - 2))Gij{4>) - h(ikGj)kict>) , (3-10) 
Hence we have an exact zero mode which may be obtained from (3.4) 

Go,ij(0) = Sij - \r]ij - Vij{4>) . (3.11) 
In consequence for (3.6) to be solvable we must have 

j d^0e-iM^^-i^-(''^)-2^(^)(5,,-i,,,-F,,(^))(r,,,-^ (3.12) 

In order to obtain tractable equations we impose 0{N) symmetry so that we need only 
deal with functions of p = and in this case we must have rjij = rjdij. Writing F((?!>) = v{p) 
(3.1) becomes 

2pv"{p) -id-2 + v)pv'ip) = -Nv'ip) + 2pv'{pf - dv{p) . (3.13) 
At the origin we must have Nv'{0) + dv{0) = and asymptotically 

v{p) l{2-r])p + Apdd=^ . (3.14) 

To ensure non singular solutions as before it is necessary to fine tune v{0). For critical 
exponents we consider, if > 1, spherical harmonics which are expressible in terms of 
symmetric traceless tensors or equivalently 

Fii(t>) = {t ■ 4>y flip) , t^ = 0, 1 = 0,1,2,.... (3.15) 

The eigenvalue equation (3.2) becomes 

I^iMp) = (Az + lUd -2 + v)- d)fi{p) , (3.16) 

for 

-Di = 2p^ -{d-2 + r/)p^ + (TV + 2/)^ _ Apv'{p)-^^ 2lv'{p) . (3.17) 
The relevant boundary conditions arc that fi{p) is analytic for p ~ and non singular for 

_d-A|__ 1^ 

p > and oc p''+^-v 2 as p ^ 00. Corresponding to (3.4) 

fiAp) = U^-ri)-v'{p), (3.18) 

satisfies 

2^i/i,<^(p) = 0, (3.19) 

so that Ai,(^ = ^{d + 2 — rj). The scalar product, with respect to which T>i is hermitian, 
becomes from (3.5) 

{fl,gi)l = / dpp2^+'-ie-2(<^-2+'?)p-Mp)/,(p)5,(^), (3.20) 
Jo 
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When = 1 only / = 0, 1 are relevant, corresponding to even,odd eigenfunctions. In terms 
of v{p), which solves (3.13), (3.12) becomes 



r/ / d/>p2^-ie-2('^-2+'?)^-2''W(l(2-r/)iV-Ar^;'(p)-2pi;"(p)) 
Jo 



— /°dpp§^-ie-^('^-2+^)^-2''(^) 
Jo 



X (i(2 - v){Nv'ipf + WipW'ip) + ^p'v"ipf) 

- Nv'{pf - 6pv'{pfv"{p) - I2p'^v'{p)v"{pf - Sp^v"{pf^ . (3.21) 

When A'' = 1 it is easy to see that this is identical with (2.19) since then V"{(p) = v'{p) + 
2pv'{p). 

To decompose (3.6) wc write 

Zij{4>) = 6ij z{p) + (t>i4>j {z'{p) + y{p)) = % {(t>j)z{p)) + MjVip) , (3.22) 

and using 

VZij{<p) - 2V^ik{<P)Zj)k{<P) = 5^{^iz{p) + 2y{p)) + 4>i(t>, ^(^^1^(P) + 2y{p)) 

+ cPicj)j {Viy{p) -{d-2 + r, + Av'{p) + Apv"{p))y{p)) , 

(3.23) 

we may reduce (3.6) to 

2d 

V,z{p) = ri- 2y{p) - — v'{pf , (3.24) 

and 

V,y{p) -{d-2 + r, + Apv"{p) + Av'{p))y{p) = pv"{p? ■ (3.25) 

The equation for y thus decouples from that for z so that (3.25) may be solved and the 
result used in (3.24) to then determine z. Asymptotically z{p) approaches a constant, just 
as in (2.10), but y(p) vanishes. Since, as shown in (3.19), Pi has a zero mode /i,,/,, given 
by (3.18), T] must be fixed to allow a non singular solution by, 

ri {h,4>, l>i = 2 </i,<^, y), + ^ </i,^, v'')^ . (3.26) 

Using 

X {Viy{p) - (d - 2 + 77 + 4.pv"{p) + 4.v' {p))y{p)) 

= -2{d-2 + ri){h,4>,y)^ (3.27) 

we may also obtain 

{d-2+r^){h,4>,y), = ^^1ippi^+ie-i('^-2+'')''-2^(^) {\(2-r^)-v\p)-2pv'\p)y {p)\ 

(3.28) 
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Combining (3.26) and (3.28) is equivalent to (3.21). If we restrict to scalar fluctuations, 
without any harmonics, then we may decompose Gij{(f)) in terms of g{p),h{p) in a similar 
fashion to (3.22) so that (3.8) reduces to 

i^ig{p) = Hp) , (3.29) 

with ^ = and also 

Vih{p) - 4{pv"{p) + v'{p))h{p) = {X + d-2 + ri)h{p) , (3.30) 

when (T>i — X)g{p) = —2h{p). Manifestly the eigensolutions for g in (3.29) are identical 
with the solutions of (3.16) for I = 1. The eigenvalues are related as expected from (1.1) so 
the solutions of (3.29) represent redundant operators. The eigenvalues determined by the 
/i-equation (3.30) give exponents corresponding to genuine physical scaling operators. 

Numerically (3.13) and (3.16) can be solved straightforwardly for d = 3, as before after 
precisely tuning /c^^^ = ?;(0) to ensure non singular solutions for all p > 0. The eigenvalues 
are denoted by where Aq^^ = 3 and we take m = 0, 1,2, . . . . In terms of the single 

component results Aq"^,^^ = A2m and A^^^ = \2m+i- When / = some results are given in 
table 3, for = 1 they match as expected the results given earlier. For r/ = these agree 
with results in [27, 28]. In the large N limit, taking = 0, then A^'^^ = 3 — 2m — 



rj 






A;(2) 


^0.1 




a(^) 

^0,1 


fc(lO) 


x{10) 
^0.1 


0.00 


0.076200 


1.5395 


0.19934 


1.4120 


0.53876 


1.2432 


1.7637 


1.0886 


0.01 


0.073513 


1.5409 


0.19176 


1.4175 


0.51693 


1.2535 


1.6929 


1.0999 


0.02 


0.070886 


1.5421 


0.18439 


1.4229 


0.49578 


1.2634 


1.6242 


1.1112 


0.03 


0.068319 


1.5433 


0.17723 


1.4281 


0.47534 


1.2732 


1.5577 


1.1224 


0.04 


0.065811 


1.5444 


0.17026 


1.4330 


0.45552 


1.2823 


1.4931 


1.1336 


0.05 


0.063360 


1.5454 


0.16349 


1.4378 


0.43633 


1.2921 


1.4306 


1.1447 



Table 3: First / = Eigenvalues for AT = 1, 2, 4, 10 



V 




x(2) 
^0,2 


^0,2 


x(10) 
^^0,2 




X(2) 
^^0,3 


^^0,3 




0.00 


-0.6557 


-0.6712 


-0.7338 


-0.8713 


-3.1800 


-3.0714 


-2.9400 


-2.8985 


0.01 


-0.6449 


-0.6596 


-0.7187 


-0.8515 


-3.1560 


-3.0498 


-2.9199 


-2.8728 


0.02 


-0.6341 


-0.6480 


-0.7039 


-0.8319 


-3.1319 


-3.0281 


-2.8998 


-2.8473 


0.03 


-0.6232 


-0.6364 


-0.6892 


-0.8124 


-3.1076 


-3.0062 


-2.8795 


-2.8221 


0.04 


-0.6123 


-0.6247 


-0.6746 


-0.7930 


-3.0832 


-2.9841 


-2.8592 


-2.7970 


0.05 


-0.6013 


-0.6131 


-0.6602 


-0.7739 


-3.0585 


-2.9617 


-2.8397 


-2.7722 



Table 4: Second and Third I = Eigenvalues for iV = 1, 2, 4, 10 



For I = 1 then A^ q = ^{5 — Vj) and A^ ^ = \{l + vi) which is a useful check. Some other 
results are given in table 5. We also present some results for Z = 2 in table 6 and I = 4 

in table 7. An important observation is that A4^g < whereas A^^'' > for N > ?>. This 
reflects the instability of the 0{N) symmetric fixed point against RG flow to one with just 
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v 


^1,2 


x(2) 
^1,2 


^1,2 


x(io) 

'^1,2 


^1,3 


x(2) 
^1,3 


a(^) 

/^l 3 


x(io) 


0.00 


-1.8867 


-1.7986 


-1.6741 


-1.5535 


-4.5244 


-4.3251 


-4.0185 


-3.6719 


0.01 


-1.8696 


-1.7835 


-1.6615 


-1.5406 


-4.4932 


-4.2984 


-3.9983 


-3.6529 


0.02 


-1.8524 


-1.7683 


-1.6487 


-1.5277 


-4.4618 


-4.2714 


-3.9778 


-3.6339 


0.03 


-1.8351 


-1.7530 


-1.6358 


-1.5148 


-4.4301 


-4.2442 


-3.9570 


-3.6149 


0.04 


-1.8177 


-1.7375 


-1.6228 


-1.5020 


-4.3982 


-4.2166 


-3.9358 


-3.5958 


0.05 


-1.8002 


-1.7220 


-1.6096 


-1.4891 


-4.3661 


-4.1887 


-3.9143 


-3.5767 



Table 5: / = 1 Eigenvalues for AT = 1, 2, 4, 10 



V 


X(2) 


^2,0 


x(io) 


X(2) 
^^2,1 


^2,1 


x(io) 

^^2,1 


0.00 


1.7819 


1.8476 


1.9283 


-0.4737 


-0.3332 


-0.1531 


0.01 


1.7760 


1.8396 


1.9187 


-0.4675 


-0.3308 


-0.1535 


0.02 


1.7701 


1.8317 


1.9091 


-0.4612 


-0.3282 


-0.1539 


0.03 


1.7642 


1.8239 


1.8996 


-0.4548 


-0.3255 


-0.1541 


0.04 


1.7583 


1.8161 


1.8901 


-0.4482 


-0.3226 


-0.1542 


0.05 


1.7525 


1.8084 


1.8807 


-0.4416 


-0.3195 


-0.1543 



Table 6:^ = 2 Eigenvalues for A?" = 2, 4, 10 



V 


X(2) 


X(3) 
^4,0 


^4,0 


^^4,0 


X(2) 
^4,1 


x(3) 
^4,1 


^4,1 


-^4.1 


0.00 


-0.0337 


0.1109 


0.2315 


0.6045 


-2.6147 


-2.3954 


-2.2093 


-1.6169 


0.01 


-0.0358 


0.1046 


0.2218 


0.5871 


-2.6023 


-2.3889 


-2.2077 


-1.6286 


0.02 


-0.0377 


0.0985 


0.2124 


0.5700 


-2.5897 


-2.3821 


-2.2056 


-1.6373 


0.03 


-0.0395 


0.0926 


0.2032 


0.5530 


-2.5768 


-2.3749 


-2.2031 


-1.6456 


0.04 


-0.0412 


0.0869 


0.1943 


0.5326 


-2.5636 


-2.3673 


-2.2001 


-1.6536 


0.05 


-0.0428 


0.0815 


0.1856 


0.5197 


-2.5501 


-2.3593 


-2.1967 


-1.6612 



Table 7: Z = 4 Eigenvalues for AT = 2, 3, 4, 10 



discrete cubic symmetry when > 3. That the critical Nc < 3 has been shown in very 
detailed multi-loop Feynman diagram calculations [29]. 

We may also use (3.21) and also (3.26) with (3.28) to determine r/^^^ when d = 3 with 
the results 

= 0.0413 , 77(2) = 0.0414 , r/(^) = 0.0390 , r/^^^ = 0.0357 , 
77^^°) = 0.0200 , r/(^°) = 0.0125 , (3.31) 

falling off as expected with increasing N. 

In order to verify consistency with large N results we also calculated eigenvalues for 
AT = 20 and r? = 0, when k^^^^ = 3.8727448, which are given in table 8. For I, m small these 
are closed to the asymptotic values. 

An interesting special case is N = —2, which was considered in [16]. For this N from 
(3.13) we get v'{0){v'{0) — 1 + ^r/) = 0. When v'{0) = 1 — ^77 we have the high temperature 
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m = 





1 


2 


3 


4 


5 


,(20) 


3.000 


1.041 


-0.937 


-2.938 


-4.966 


-7.025 


,(20) 


2.500 


0.500 


-1.521 


-3.566 


-5.640 


-7.744 


J20) 
^2,m 


1.963 


-0.077 


-2.139 


-4.227 






,(20) 
^3,m 


1.392 


-0.687 


-2.789 


-4.917 






,(20) 
^4,m 


0.788 


-1.328 


-3.468 


-5.634 







Table 8: Eigenvalues for AT = 20 



fixed point solution v{p) = (1 — |r/)(p+ |). The relevant non trivial fixed point arises for 
v{0) = v'{0) = 0. Consistency with our equations requires r]^~^^ = 0. To show this we may 
note in (3.18) that /i,</,(p) — >■ ^(2 — 77) as p — so that, with the definition of the scalar 
product in (3.20), the integral (/i,</,,l)i diverges due to the singular behaviour as p — > 0. 
On the other hand and also given by (3.28), are finite. Hence (3.26) is 

only consistent when r/ = 0. For the corresponding eigenfunctions it is necessary to consider 
boundary conditions carefully. In (3.16) then as p ^ the differential equation requires 
solutions for fi{p) which are 0(1) and 0{p°-) where a = 1 - ^{N + 21) - 2Nk/d, k = v(0). 
Generally we impose the requirement that the second solution is absent. For A'^ = —2, I = 0, 
when also k = 0, there are solutions with fo{p) = O(p^) for which the norm given by (3.20) 
is finite. However to obtain eigenvalues which are related to those for A/" > it appears 
necessary to consider as before solutions with /o(0) = 1, and we may impose /o"(0) = 0. In 
this case we find numerically Ao,i = 2. 



4 £ Expansion 



A further consistency check, following [21], is to consider the £-expansion where 
2n d 



n — \ 

In this case we define 



d-2 



n+Un-lf£ + 0{£^), n = 2,3,.. 



x = Ud-2)p, 



and, with v{p) = v{x), rj = 0, (3.13) becomes 

X v"{x) + {a + l-x) v'{x) + — ^ v{x) = X v'{xf , a=\N -1. 

Using (4.1) it is easy to see that as £ — there is a solution of the form 

v{x) = kneK{x) + 0{e^), 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



for a Laguerre polynomial, satisfying xV^' ix) + (a + 1 — x)L'^{x) + nL"(x) = with 
the orthogonality condition 



/ dx e-^x" V^{x)Vl{x) = ^ r(a + 1 + r) 6rs ■ 
Jo f- 



(4.5) 
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To determine kn we consider the scalar product of both sides of (4.3) with to 0(£^) 
giving 

rco f'CO 

l{n-lfkn dxe-^x"L^{xf = kr? dx e-^x"+i L^(a;) L^'(a;)2 . (4.6) 
Jo Jo 

This requires 

^,(_.)„ (n-l)^(« + l)„ 

nnl ry(<^:^) 

where we define 

= v( Li^u / ^?^'(^) ^^'(^) ^*"(^) ■ (^•^) 

i (a + / + Ij Jo 

It is straightforward to set up a perturbation expansion in e for the higher order terms in 
the solution (4.4) as a series summing over Laguerre polynomials L°. 

For the associated eigenfunctions and eigenvalues then letting 

\rn = d-{d-2){m+\l + \rn), (4-9) 

from (3.16) and (3.17) we require that Xi^m is determined by 

Vlfl,m{x) = - {km + m) fl,m{x) , (4.10) 

where, with 77 = and a as in (4.3), T>i = {d — 2)'Di, 

d^ d d 

Vi = x—^ + (a + l + l-x)- 2xv'(x)- Iv'(x) . (4.11) 

ax'^ ax ax 

Note that we must have Ai,o = 0) M,n-i = ^{n ~ — 2) to ensure the exact results 
Ai,o = \{d + 2), \\^n-i = \,{d — 2). As £ — it is easy to see that 

fi^m{x) ^ L^+'(x) , A,,^ = 0(£) . (4.12) 
To first order, where we may use (4.4) for v in (4.11), perturbation theory gives 

j^dx e--x°+' L-^\x) {2xL-\x)^ + I L«'(x)) L'^\x) 



I _ , ml 



-(_-i\n , Wt! Cmmn . s 2 (« + l)n ml G^^;^n , . 

"^-(a + ^ + l)^^-(" ^K\G^J^^{o. + l + lU'- 
To lowest order r/ = 0{£^). Using (3.26) and (3.28) with (4.1) and (4.4) gives 

poo n r, 2 2 / /-oo /-oo \ 

ri j dx e-^x"+i = -J^^Zr^ ( y e-^x"+iL^'(a;)2 + 2 j dx e-''x''+''L'^^'{xf\ . 

(4.14) 

Noting that L'^'{x) = —L'^'^\{x) we then find from (4.5) 
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Explicit results in (4.13) and (4.15) can be obtained using 

^(a,l) ^ {a + l + l)r {a + l + l)s 
~ {l + r + s-ty. 

^ V- 1 {l + 2n)l 1 ,^ 

^ n\{r-n)\{s-n)l {t-r-s + 2n)\ {a + l + l)n' 

Prom (4.16) we have G^i'^ = Q so that G^i^^ = 0, t > 2, and (a+l)G£:i^„_i„ = 

These ensure that Ai^o = 0, Xi,n-i = ^{n — 1)^£ as required. For the case of primary 

relevance here n = 2 and we have 

A/,m = ^2 {{I + 2m){l + 2m - 1) + m{N + 21 + 2m - 2)) e . (4.17) 

These results are in accord with those obtained [30, 31] in early calculations involving the 
e-expansion. For the potentially marginal operators A4^o = 2(N — 4)k2£, \2,i = — 14fe2£ and 
Ao,2 = so that the critical A'^c = 4 to leading order. 

We may also extend the 0(e) calculations to (3.30) for the non redundant derivative 
operators. With t] = and the same change of variables as in (4.2) this becomes 

(Pi - 2xv"{x) - 2v'{x))hm{x) = - {m + Xm)hm{x) , 

Am = -(d-2)(m + l + A^), m = 0,l,..., (4.18) 

where with (4.4) we may take hm{x) = L^+^(x) + 0(e) and A^ = 0(£). To first order 
perturbation theory gives 

Xm = Xi,m + 2Ke • / dx e--x"+' L^+\x)^xLr{x) + L^'ix)) 

= Ai,„ + {-ir2k^e ((n + a) gL";^Li - « G^^^i^^^) , (4.19) 

using the identity xL^"{x) + L^'{x) = — (n + a) L"_^(a;) + aL^'tl{x). When n = 2 this 
gives 

Am = Ai,m + k2e{N + 8m + 2) , Ai,™ = m{N + 6m + 2) , (4.20) 

giving 

J (m + l)iV + 6m2 + 10m + 2 



5 Perturbative Discussion 

The crucial significance of the e-expansion is that it is possible to use standard perturbative 
methods in quantum field theory involving a Feynman graph expansion. To show the 
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parallel with the above treatment we describe how the same results are found in terms of 
conventional /3-functions. For a general scalar lagrangian 

i5'^0-5^<^ + y(<^), (5.1) 

then we may define a generalised /3-function for the potential V, which is a linear function 
of the couplings^ , where 

Bvicf>) = - 2)<^ • dV{cl>) - dV{<i)) + I3v{ct>) ■ (5.2) 

The perturbative RG flow equations are then 

V{(t>) = -Bvi4>) , (5.3) 

and the fixed points V = V^ are hence defined by 

BvM = 0. (5.4) 



With d as in (4.1) then for a renormalisable theory V{(f)) is a polynomial of degree 
2n. Using background field methods and minimal subtraction [32] then in a perturbative 
expansion calculations give /3y in the form 

/3y(<^) = (700) • dVicj)) + Pv{^) , (5.5) 

where ^y^^ij is the anomalous dimension matrix for the field (f) and Pvi4') depends on 4> solely 
in terms of scalar contractions of various products of ^ij2...ifc (</>) = di^ . . . di^V{(j)) with 
k > 2. For renormalisable V, 7^ depends only on the ^-independent giii2...i2n — ^ii2---»2n 
and /3v{4>) is also just a polynomial of degree 2n. In general there are contributions to /3y at 
(p- l)(n- 1) loops, p = 2,3 . . . , when pv = 0(52"(p-i)FP) and, if p > 2, -f^{g) = 0{gP-^). 

Assuming K(0) = j^g*iii2...i2„(t>ii ■ ■ ■ 4>i2r, t^ien (5.4) determines, for n = 2,3,..., 
g*hi2...i2n perturbatively as an expansion in e. At the fixed point the anomalous dimension 
matrix for cj) also defines 

In the vicinity of a fixed point defining, for -F(0) an arbitrary polynomial of degree 2n, 
a linear operator 7 by 

/?y+F(0) = (5v{^) + ^F{4>) + 0(^2) ^ (5 7) 

then critical exponents are determined by the eigenvalue equation 

AF(<^) = \{d - 2)<p . dFi<p) + 7,F(<^) = id- A)F(0) , 7* = l\v=v. ■ (^-S) 

As in the discussion of the exact RG equations there are explicit eigenfunctions relating to 
(j) itself. Corresponding to (5.5) we have 

7 = (7,^<^)-5 + 7, (5.9) 

^For V{cj)) = '^iq'Oi, where form a basis of monomials in 0, then (3v{<f) = (^^ {9)^1- Subse- 
quently/?v^=E//3'(s)a?T. 
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where 7 involves at least second order ^derivatives. Hence 

A(j)i = l{d - 2)(t)i + lriij(j)j , (5.10) 

and also by differentiating (5.4) 

^V^A'P) = Ud + 2)K,i(<^) - i%Kj(</)) . (5.11) 

C is arbitrary up to total derivatives so in this discussion mixing with operators which are 
just spatial derivatives is neglected. The relevant matrix is triangular so the determination of 
scale dimensions is not affected. Up to a derivative operator proportional to K,i(0) ~ 
0, as a consequence of the field equations so this is redundant. 

At lowest order, as shown using background field techniques with n — 1 loop Feynman 
diagrams in [21], 

/?^(<^)("-i) = anVi,...ijmi-iM) , (5.12) 
where higher loops are 0{V^). Also, for 2(n — 1) loops. 



where 



(n!)2 ^ 

if ' ~ ^ (2n)! ^"l*2--i2n-l 5'jili2--j2n-l 5 (5.13) 



an = ^^r(^r\ (5.14) 
" {AttY' n! Vn-1/ ^ ^ 

The definition (5.7) and (5.12), then determines to lowest order 

7^"-^) = 2anVi,...iM)Fh-iA<P) ■ (5.15) 

As before we impose 0{N) symmetry so that 

V{^)=g^{\<j>''r, /3y(0) = /?,(5)l(i</.2)", 7^,i,=70(5)%. (5.16) 
At the fixed point (5.4) 

(ig{g*) = (n - 1)5* £, r} = 270(5*) . (5.17) 
With 0{N) symmetry the eigenfunctions in (5.8) have the form 

Fi,M = {t-cpyc^^^r, t' = o. (5.18) 

The corresponding eigenvalues are then 

d - Xi,m = Ud-2 + r]){l + 2m) + 7,,^ , (5.19) 
with ^i^fn the anomalous dimension determined by %Fi^m{(f)) = 7z,m-^/,m(0). 

To handle the combinatorics involved in evaluating (5.15) with F as in (5.18) we re- 
express this using 

^■^) (5.20) 
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and then follow the method described in [31]. First we note 



and then 



= S rUk-Jl<t-k+rn (i-t-'r-'-^' ■ a'r"'(is"r({t ■ mi^p'^r) ■ (5.22) 



Using 



ilc/'nit ■ cj^'Yiy'r) = — ^ ^a + l + l + m-r)r{t- <t>')\yr~' , (5.23) 

[m r). 



with a as in (4.3), and 



we then obtain^ 

s (I ■ ^)'((^*')" <* ■ L.., = -^-i ■ *)'a*^)"*"-' . 

^-'- = ^r!(fc-2r)!(n-A: + r)! (m - r)! + ^ + + " (2m + Z-A;)! ' 

Hence from (5.12) and (5.15) with (5.16) and (5.18) we may obtain to lowest order 

/3,(5)("-^)=4n'Sa„5', (5.26) 
Prom the fixed point equation (5.17), with 77 = and (5.26), we get 

4n',°2«n5* = (n-l)£. (5.27) 

In consequence 

7i""'^i^«,m(</') = 2a„ 5* At% Fi,^{cP) . (5.28) 
and, in (5.8) and (5.19), to lowest order 

7;,^ = 2an 9* 4"ni = 2(n - 1) ^ e . (5.29) 



^The coefficients A^™'^^ satisfy various identities, in particular 

(a + I) {Atf,^^^ + 2m(a + I + m) ^i^'^i.,) 

= m (2a + 0(« + / + m + n - fe) A^-^i^'+i^ + Z (a + / + m)(m + n - fc) A<"„;-'' • 

WW AT - 1 yd^-S''^ - (2n)!(2m+0! for 7 - 1 
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Note that 7o,n = Pg{9*)- 

As a special case of (5.25) we also have 

so that (5.6) and (5.13) give 

ri = 2n! (a + 2)„_i (ung^f = 2in-lfn\ (5.31) 

The precise identity of the results (5.29) and (5.31) with those obtained from the RG 
equation (4.13) and (4.15), where A^,^ = ^(n — l)7;,m) follows from 

When n = 2 the 0(e) results may be read off from 
'^^^3* = j;rp8' 42;m = (2"i + 0(2m + l-l) + m(iV + 2/ + 2TO-2). (5.33) 

6 Mixing Effects 

For operators which are monomials 0™ with m > 2n then perturbatively it is necessary 
to include mixing effects with operators (9^)^^"*"^*^. For m > 4n — 2 there is additional 
mixing with operators involving four derivatives, such as (9^0)^, but this is neglected here. 
We here discuss how the treatment of the previous section may be extended and show how 
reparameterisation invariance is manifest in a perturbative approach. 

The initial renormalisable lagrangian £ in (5.1) is extended to 

C = Cv + Cf,g, Cv = ld^cf> -8^,(1) + V{(I)), Cf,g = F{(I,) + \G^,{<t))d>'^id^<t), . (6.1) 

Although for general F, Gij the resulting C is non renormalisable, keeping only counterterms 
which are linear in F, Gij , we may consistently define a bare lagrangian Cq which extends 
the renormalisable theory defined by Cy to include first order perturbations by finite two 
derivative operators, so long as F{(j)),Gij{(f)) are constrained to avoid the necessity of four 
derivative counterterms. As usual there are corresponding /3-functions 

BF{<t>) = lid - 2)<t> ■ dF{<t>) - dF{<t>) + PF{<t>) , 
BG,iM) = \{d - 2)0 • dGij{(t>) + i4>,ikGkj{(t>) + l<t>,3kGik{(t>) + fe(0) , (6.2) 

which are linear in F, Gij so that 

f3p = JFPF + jFG,ijGij , PG,ij = lGF,ijF + jGG,ijklGkl ■ (6.3) 
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Here 'yFF:7FG,ijilGF,iji7FG,ijkl are differential operators depending on the renormalisable 
couplings or V, clearly we have, restricted to F{(p) of degree 2n, 'ypp = 7 as defined in 
(5.7). At a fixed point the exponents are defined by the coupled equations 

BF{<l>)\y^y^ = -\F{<P) , BG,ij{cP)\y^y^ = -XGiM) ■ (6-4) 

For the lagrangians in (6.1) we have 

6Cv = dF,G^F,G for 6(l)i = Vi{(j)) , 6dn(l)i = Vij{(f))d^(l)j , (6.5) 

if3 

SF,GFi(t>) = vi4>)-dV {(!)), SF,GGiji4>) = diVji(t>) + djViict)) . (6.6) 

If F{(p), Gij{(f)) are restricted to ensure that no mixing with four derivative operators arises 
then it is necessary to require Vi{(f)) = O(0^("~^)). 

As was apparent in the discussion of renormalisation for general two dimensional a- 
models [33, 34] invariance under reparameterisations S(pi = Vi{(j)) leads to a corresponding 
freedom in the definition of the /3-functions. Here we show how this leads to relations for 
the exponents defined by (6.4). Assuming first 

f^F{<P)\F=..av,G.,=a.v,+d,.. = ■ dV{<P) + v{^) . dPyicP) , (6.7a) 

PG,ij{<P)\F=vdV,GM=dkvi+diVk = - '^l4>,jkVk{^)) + djilvM) - '^l<l>,ikVk{(t>)) , 

(6.7b) 

with 7 defined in (5.8), then (6.2) gives 

BF(^\F=..av,G.,=a..,^a,n = ^^^^ ' ^^^^^ + "('^) " ^^^('^) ' (^.Sa) 

BG,M)\F=..av,G,^=a,n^d^v, = ^^^M) + d^m) , (6.8b) 

for 

[/,((/.) = i(d - 2)0 • dViicP) + -fV^{(P) -lid- 2)V^{^) - -f^,^JVj{(P) . (6.9) 

A general justification of (6. 8a), (6. 8b) with (6.9) is described in appendix C. 

At a critical point, where (5.4) holds, then for vector solutions of (5.8), Auj = {d—\y)vi, 
there are, as a consequence of (6. 8a), (6. 8b) corresponding solutions of (6.4) such that 

\ = K-\{d + 2-ri), (6.10) 

assuming a diagonal form for r]ij. Thus in this perturbative context we reproduce (1.1). In 
particular as a consequence of (5.10) we have the exact zero modes 

Fo(</.) = • dV{(t)) , Go,,,- = 25ij . (6.11) 

^The full Lagrangian in (6.1) is invariant if (6.6) is extended to 6F = v ■ d{V + F), 6Gij = diVj + djVi + 
V ■ dGij + diVk Gkj + djVk Gik- 
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To verify these results we consider the lowest order perturbative results at n — 1 loops 

^ — ^ 77-! 

r,s,t>l 
r+s+f— n+1 



r>2,s,t>l 
r-\-s-\-t=n-\-l 



^ — 77,! 

Z-^ r\ s\t\ ^h-irki...kti<P)^ji...jski...kti'l^) <^ni2,i3-vji-is (^) > 



(6.12) 



where a„ is as in (5.14) and 



r(^^)r(^)r(^) 
^"-^^r(^)r(^)r(^) ^ ^ 

When r = l,s + t = n, Krst = 1- Also 

- ^jl...j„(j(0) G'iji2j)i3- • (6.14) 

We may directly verify that (6.12) and (6.14) satisfy (6. 7a), (6. 7b) with -y^"^-"^) given by 
(5.15) and also 7^ = 0. At the order given in (6.14) there are no contributions involving F. 
At the next non zero order 

+ (2n - l)4a„2 |^ gi,...i,r.-2kii 9^,...^2n-2lJ) Gkl{(p) + 0{dG) , (6.15) 

which is also compatible with (6. 7a), (6. 7b) using (5.13) for 7^. 

For definite results we assume 0{N) symmetry as in (5.16). To first order in e the 
equations (6.4) decouple and we may write the eigenvalue equation 

i(d - 2) ■ dGij{(l)) + 2an9* 'DGijicj)) = -A (<^) , (6.16) 

where from (6.14) T)Gij is given by 



n\ \d(j) d(j)' 



2/(9 d 

4>'=4, n! 



Clearly for Gij{(f)) = Sij then A = 0. In general we have 

- Ju^i^^^il^'^'"') ^ (6.18) 
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so that from (5.25) 

GiM = Sijil4>T+^ + (m + 1) (l>i(l>jil4>T , m = 0, 1, . . . , 

A„ = - (d - 2)(m + 1) - 2a„5*4n,m+i = '(d - 2)(m + 1) - 71,^+1 , (6.19) 

as expected from the general relation (6.10) with (5.19) and 77 = 0. In general 
1 



^ — V{S,,i'^c{>^r+') = am Sij{l<PT+' + Pm Mj{l<f>T , (6.20) 



and we then find another eigenvalue, in addition to (6.19), which can ne expressed in the 
form 

Am = -(d - 2)(m + 1) - 7m , 7m = 2an5*((m+ l)am -/5m) • (6-21) 
To obtain more explicit results for this we apply (6.17) to obtain 

+ 2((m + l)^;r_:?„,m - 2mAl":l)„,„,_i - (« + n)m4!:^i-i,m-i) 

X (,5,,(i02)m^^^.^.(1^2^m-l) 

1 (JL . JLy~)(1 



+ + (^r^ (I ■ |^)""«^*''""' ^'"'■^y^r) . (6.22) 
Prom this we may determine 

(m + l)am - /3m = ^l"',m+l 

+ ^ (m + l)(a + m + 1) {A^:^^^^ + (a + n)A^::Sn-iJ 
4 



--m{a + m + 2) (4-iU-i + (« + ^)4-l-i,m-i) • (6-23) 

Combining (6.23) with (6.21) gives the exponents for new non redundant operators so long 
as > 2, for N = 1 the corresponding eigenfunction vanishes. When n = 2 we have 

(m + l)am - I3m = {m + 2)N + 6m^ + 12m + 4 . (6.24) 

Hence combining this with (5.33) the anomalous dimensions of these derivative operators 
at the fixed point, to first order in e, are given by 

(m + 2)Af + 6m2 + 12m + 4 
7m = e. (6.25) 

This may be compared with 2Am in (4.21) which was obtained from the approximate deriva- 
tive expansion for the Polchinski RG equations. Although similar they are not identical. 
The perturbative results in (6.25) are of course the first terms in a well defined expansion 
to any order in powers of e. 
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7 Conclusion 



The status of the equations presented in this paper for extending the local potential ap- 
proximation to the Polchinski exact RG equation is unclear, in that there is no consistent 
derivation and the resulting equations for V, Z are partially decoupled in that the F-equation 
lacks expected Z contributions. Nevertheless the LPA equation now involves ij which, in this 
respect, is similar to an approximation made in a treatment of exact RG equations in [35]. It 
remains to be seen whether the introduction of terms involving Z into the Polchinski LPA 
equations is at all possible, while maintaining the crucial property of reparameterisation 
invariance, and so allowing a well defined determination of 77. Some time ago Morris [16] 
obtained, with a particular cut off function K{p^) proportional to a simple power ofp^, and 
for the exact RG equations applicable to the one particle irreducible functional F, a set of 
equations, in a derivative expansion, which are invariant under global rescaling of the fields. 
To this extent reparameterisation invariance is preserved and there is consequently an exact 
zero mode so that rj may be determined unambiguously. However these equations are highly 
nonlinear and hard to analyse. The associated zero mode eigenfunction for these equations 
has not apparently been found explicitly in the literature. If a derivative expansion is to be 
consistent then results should not change dramatically when the LPA for V is extended to 
first order in derivative operators to a pair of coupled equations for V, Z. However results 
for some eigenvalues obtained using the Morris equations differ significantly [26]. It would 
be very desirable to understand more analytically what features of the equations in [16] 
ensure reparameterisation invariance so that this crucial constraint might be imposed more 
generally in derivative expansion RG equations. It might also be helpful, as was the case 
here, to construct an appropriate scalar product for the eigenfunctions whose corresponding 
eigenvalues are the essential output of exact RG equations. 

In general reparameterisation invariance, which is related to issues of scheme indepen- 
dence, is akin to a gauge invariance of the exact RG equations [36]. As is well known for 
gauge theories violating gauge invariance in some approximation can lead to unphysical 
conclusions. Emphasising the consequences of reparameterisation invariance, and the con- 
sequential presence of an exact zero mode in the RG flow equations, might also be a useful 
criteria for optimisation in equations where reparameterisation invariance is not automatic. 

It is of course non trivial that the LPA captures the global aspects of RG flows between 
various possible fixed points in scalar field theories. The calculated critical exponents are 
also of essentially the correct magnitude since the results must agree with those from ordi- 
nary perturbation theory for non derivative operators to first order in e (strictly this appears 
to have been shown only when maximal 0{N) symmetry is required, it is presumably true 
for fixed points with lesser symmetry although a general argument appears to be lacking) . A 
natural constraint for any derivative expansion is that this agreement should extend to non 
redundant scalar operators with two derivatives. The discrepancy between the anomalous 
dimensions for such operators given by (4.21) and (6.25) shows that the equations discussed 
here are not fully satisfactory in that respect. Perhaps more general derivative expansion 
equations can be obtained by using more input from perturbative results. 
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A Large N Calculations 



Solutions of the LPA fixed point equations become more tractable in the large limit 
[27, 37]. We show here how these results may be used to obtain a corresponding value for rj 
using (3.21) or (3.26) and (3.28). In the large N limit we suppose v,p = 0{N) and taking 
the derivative of (3.13) then gives 

{{d-2 + ri)p- N + Apv')v" = (2 - r])v' - 2v''^ , (A.l) 

or 

2v'{s-v')^-{d-2s + iv')p = -N, s = l-ir/. (A.2) 
Such first order linear equations are readily solved giving /> as a function of v', 



with the integral extended by analytic continuation for > 1. To ensure a smooth 
continuously differcntiable solution for all p > wc must set the constant of integration 
C = and then the large N result can be expressed as 

p=^^(2,i;2-^rf;y), (A.4) 

for F a standard hypergeometric function. 

In the large N limit ry — so we may set s = 1. When d = 3, F{2, 1; |; u') = for 
v' = -0.6349132 and so ^(0) ~ 0.2116377A/". For general d we may invert (A.4) giving 

v'ip) = iP - Po) + 0((p - po?) , Po = ^ . (A.5) 

In integrals obtained from the scalar product (3.20) we then have 

p\N ^-Ud-2)p-2v{p) ^ce-'2w iP-Po? , (A.6) 

for some constant C, and so the dominant contribution for large N arises for p ^ po where 
we may use (A.5). Hence 

{fi,^,v"), ^ i4-d)\d-2) 

Furthermore in (3.28) for p po 

p(i (2 - r;) - v'ip) - 2pv"{p)) « \N , (A.8) 

so that (3.28) gives 

(A.0,y)i ^ d{d-2){4-df 

</i,0,l>i ^ 8{d + 2)N ■ ^ ■ ^ 

Hence from (3.26) the leading large N result for rj in this analysis is determined to be 

The exact leading order large N result 77 ~ 2(4 - d)r{d - 1) sin7r5(d - 2)/{'n:dr{ldfN) is 
numerically close to (A. 10), coinciding as d — 4. 
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B Perturbative Calculations 



We here outline how the perturbative results (6.12) and (6.14) are obtained following the 
background field methods, with a background field <^ and dimensional regularisation, used 
in [21]. To obtain (6.14) we consider vacuum graphs which are 0{VG) 

W"i = X]^ / d'^a;id'^a;2 {Gij,n...v(¥?i)9''¥?ii9^</'ijGo(a;i2)'' 

+ r(r - 1) Gi,i,,i,...iM 5''Go(xi2)5^Go(xi2) GoixuY-^]Vi,...iM , (B.l) 

for (fii = ipi{xi) etc., and Go{x) the basic d-dimensional propagator —d^Go{x) = S'^{x). 
Using, with d as in (4.1) and a„ as in (5.14), 

G,(xr - - ^ S\x) , d,Go{x) GoixT-' ~ - ili^lll^ d^5<^(a:) , 
e n — 1 en — 1 

d'^Go{x)d^Go{x) Go(x)"-2 ^ ^ 52^d(^^ _ (B_2) 

e n — 1 

Consequently the necessary counterterm to cancel the e-pole in (B.l) is given by 

^ic.t.(¥') = ^ d'^(pidu,(pj{Vi,.„iM Gij,n...in{'P) + 2T^n...i„ »(</?) Gji,^i^.„i^{(p) 

M) , (B.3) 

which directly gives (6.14). 

Correspondingly to obtain (6.12) we consider vacuum graphs which are 0{V'^G) 

W2 = -2YI ^:[^_ J d'^xid'^X2d'^X3Gij^i^,,s,j,...j,{ipi)Vii,,,,i^ki...kt{V'2)Vj^^ 
T>-i° X di'Goixu) GoixuY d^Goixu) Goix^Y Go(x23)* 

s% xdt'Go{xi2)d^Go{xi2)Go{xi2YGo{xi3YGo{x23Y- (B.4) 



By considering the pole in {X12) ^^{xis) ^^(^23^) when Ai + A2 + A3 = d, assuming 
Ai, A2, A3 < ^d so that there are no sub-divergencies, we find, for Kj-st as in (6.13), 



a^Go(xi2) Go(xi2)'' d^Go{xi3) Go{xisY G'o(x23)* 



2 nlttn Kr+ls+lt \rd 



£ n-l (r + l)(s + 1) 
5'^Go(xi2)5^Go(xi2) Go{xi2Y GoixuY ^0(^23)* 



r+s+t=n—l 

S^xu)S''{x,s) , 



r+s+t=n—l 

2 n!a„ Kr+2st ^d/ \rd/ \ m K^ 

T ^^V^vTTV^ ^ (^12)') (a;i3) ■ (B.5) 

£ n — 1 (r + 2)(r + 1) 
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Hence (B.4) requires the counterterm 



^ r,s,t>l 
r-\-s-\-t=n-\-\ 

X] ^:|Tj7| -^'•s*^n---irfei...fet(</^)^ji...ji<fe^ ' 

T->2,s,t>l • • • / 
r+s+t=n+l 

(B.6) 

which leads to (6.12). 

A consistency check may be obtained by considering (B.6) when Gij is a constant. The 
result is then C2c.t.{'f) = ~\ ^iii..i„-i (v')^ii...i„-i (v')G'.y . This may also be obtained 
directly from lowest order counterterm for V using that the kinetic term in C defines a 
metric 5ij + Gij. A similar argument is sufficient to obtain the last line of (6.15). 



C Analysis with Dimensional Regularisation 

Using dimensional regularisation with minimal subtraction for a theory defined by (5.1) and 
(6.1) the diffeomorphism invariance under (6.6) extends to the regularised theory and we 

demonstrate here the consequences for the perturbative /3-functions. The bare lagrangian, 
including all counterterms involving poles in e necessary for finiteness to first order in F, Gij , 
has the general form 

Co = f^-'{hd>'<t> ■ Zd^<t> + V(</)) + H^) + Qij{<t>) \d>'<t>id^ct>3) > (C-1) 

where ^ is the regularisation mass scale and J^, Qij are linear in F, G^i and Zij = Zji depends 
only on gi^...i2„- The bare couplings and field are defined so as to absorb all dependence 
on II. For the standard renormalisable theory, when !F, Qij are zero and V((/>) is just a 
polynomial of degree 2n, the /3-functions and 0-anomalous dimension are defined through 

(-e-(7.<^)-^ + /3v-|,)v(^) = 0, 

- e + f3v ■ ^ j Zij = '-f<p^ikZkj + ^<j,,jkZki , (C.2) 
where in terms of the /3-functions considered earlier 

M<P) = -l£(l)-dV{cj))+eV{<P) + pv{'P), %i3 = -\e5^,+^^^i,. (C.3) 
For the extended theory, keeping only contributions to first order in F, Gij , we also have 

( - £ - ■ — + Vp)H^) = , (C.4a) 

( - £ - ■ ^ + ^Z?) = l4>,ikQkM) + l4>,jkQki{4>) , (C.4b) 
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for 

P, = 3^-^ + 3.-|. + /^G,.r J^. (C.5) 

Here Pf is related to Pf just as $v is to Pv in (C.3) and PG,ij{4') — (t^'dGij{(j))+l3G,ij{4>)- 
For transformations as in (6.6) invariance of the regularised theory requires 

^iji^)\F=vdv,Gki=dkVi+diVk = '^vGij{(l>) = diVk{(p)Zkj + djVk{(f))Zki , (C.6b) 

for 

= (i; • 5y) • — + (Sfct;; + divk) ■ ^ . (C.7) 

In (C.6a),(C.6b) Vi{(p) = Vi^cp) + . . . where we may have higher order terms involving poles 
in e. 

The crucial constraints arise from the consistency conditions between (C.4a),(C.4b) and 
(C.6a),(C.6b) giving 

[D,,V^]T{(I)) = u{^) • dV{cl)) , (C.8a) 
[Dy, Vp] Gij{(p) = diUk{(t>)Zkj + djUk{(p)Zki , (C.8b) 

where 
and 

To relate to v{^) we require that all counterterms are determined by Cq. To achieve 
this we assume 

V^{<t^) = V^■-^m. (C.ll) 

Using (C.IO) this leads to 

9 ^ . \ 9 / d \ d 

^kl 

Comparing (C.8a),(C.8b) and (C.9) with (C.6a),(C.6b) and (C.7) requires 



Vy(iF{4>) - v{(l)) ■ dPvU^) = u{(P) ■ dV{4>) , (C.13a) 

T^vPgM^) = ^^UJ{<P) + djUiicP) , (C.13b) 

for 

d - 

Uiicp) = '"i- Pf{(P) - £Vi{<i)) - %,ijVj{4>) . (C.14) 
These results are equivalent to (6. 8a), (6. 8b) with (6.9). 
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D Singularities of Solutions of RG Equation 



The critical requirement for solving the RG equation (3.1) is that it is necessary to fine tune 
k = v{0) to ensure there are no singularities for any real positive p. Nevertheless there are 
necessarily singularities elsewhere in the complex plane. As shown in [26] it is of interest 
to determine the location of such singularities so as to allow the use of conformal mapping 
techniques. The structure of the differential equation determines that the singularities are 
simple poles of the form 

^'(z^) ' (D.l) 

po e*°'^ — p 

where reality of the equation ensures that zta must both give singularities unless a = 1. 
We restrict then < a < 1. The singularities are determined by numerically integrating 
along lines of constant argument, for v{0) = k and matching with (D.l). 

When d = 3 the only singularities that are found within the radius in which the numerical 
solution is valid are on the negative real axis, a = 1. Choosing 77 = the values of k for 
suitable N are given in Table 3 and the position of the closest singularity to the origin is 
given in Table 9. 





Po 


1 


2.862 


2 


2.836 


3 


2.871 


4 


2.954 


10 


3.871 



Table 9: Values of po for AT = 1, 2, 3, 4, 10 



We have also considered d = | since there are then two solutions of physical interest 
with the additional solution arising for d < 3 and representing the tricritical fixed point. 

For the solution corresponding to the standard Wilson-Fisher fixed point the singularities 
are again just on the negative real axis. The results are in Table 10. For the tricritical 



N 


k 


po 


1 


0.252995579 


1.349 


2 


0.800566594 


1.116 


3 


1.587372474 


1.014 


4 


2.451092846 


1.002 


10 


7.725892940 


1.275 



Table 10: Values of po and A; for AT = 1, 2, 3, 4, 10 

case there are genuine complex singularities, results are given in Table 11. An illustrative 
numerical solution for v'{p) compared with the pure pole term in (D.l) is shown in Figure 
3. 
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N 


k 


Po 


a 


1 


-0.043027023 


6.3636 


0.5150 


2 


-0.132093526 


6.4736 


0.4766 


3 


-0.279532538 


6.5380 


0.4383 


4 


-0.499803930 


6.5544 


0.3993 



Table 11: Tricritical values of k, po and a for N = 1,2,3,4 
V 




Figure 3: = 1 tricritical solution for v'{ip), red and yellow lines are the real and imaginary 
parts, the green and blue lines correspond to a pure pole as in (D.l) with po from Table 11 
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